Relaxation of high-energy quasiparticle distributions: electron-electron scattering in a 

two-dimensional electron gas 
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A theory is developed for the evolution of the non-equilibrium distribution of quasiparticles when 
the scattering rate decreases due to particle collisions. We propose a "modified one-collision approx- 
imation" which is most effective for high-energy quasiparticle distributions. This method is used 
to explain novel measurements of the non-monotonic energy dependence of the signal of scattered 
electrons in a 2D system. The observed effect is related to a crossover from the ballistic to the 
hydrodynamic regime of electron flow. 
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I. INTRODUCTION 

The propagation of an electron beam in a two- 
dimensional degenerate electron gas (2DEG) can be used 
to study the role of electron-electron (e-e) collisions in 
reduced dimensions, see e.g. RefsiS^. Of main interest 
are "high-energy" beams with a typical electron energy 
that exceeds considerably the thermal energy of the equi- 
librium electrons within the 2DEG, but remains signifi- 
cantly below the Fermi energy. A special combination of 
ac- and dc-measurement techniques allows for an energy 
selective investigation of electrons traversing a defined 
distance within a 2DEG 4 -. Some of these experiments 
were analyzed in terms of the r-approximatior. 5,6 , where 
it is supposed that each collision scatters an electron out 
of the beam. However, this approach is obviously insuf- 
ficient for the case of electron-electron (ee) scattering in 
a 2DEG, where small angle scattering events dominate, 
and scattered electrons may remain within the region of 
the electron beanALSiiLiii. 

For an adequate theoretical description of ee-scattering 
it is necessary to take into account the incoming terms in 
the collision operator J (see the second term on the right 
of the equation for J) of the dynamic integro-diffcrcntial 
linearized Boltzmann equation: 



Jf = -vf 



J dpV pP '/(p')> v = J dp'v p , p , 
dt dt dr dp 



(1) 



Here, / is the distribution function of non-equilibrium 
electrons, r determines the location of a non-equilibrium 
electron, and v and p are the electron velocity and mo- 
mentum, respectively (quasi-classical approximation) . F 
is the Lorentz force acting on the electron. The function 
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FIG. 1: Scheme of the e-e spectroscopy experiment (see text). 



fp/p determines the probability of non-equilibrium elec- 
trons (t'p'p > 0) or holes < 0) to be found in state 
p after scattering from state p'. A hole in this case de- 
scribes the absence of an equilibrium electron below the 
Fermi surface. 

In the general case, the solution of the kinetic Boltz- 
mann equation is a very difficult task, both in analyti- 
cal and numerical calculations. But it is well known from 
the kinetics of gases that two regimes exist which allow 
for an essential simplification of the description. First, in 
the few-collision case the probability is small for an elec- 
tron to be scattered over a characteristic length scale L ( 
e.g., the distance between the injector to the detector, see 
Figffl below) , which is much smaller than the electron- 
electron mean-free path, I, i.e. I = vv~ l L. There- 
fore, using perturbation theory it is possible to solve the 
Boltzmann equation by iteration of the collision opera- 
tor (quasiballistic approximation). Second, in the case 
of frequent collisions (I L) a hydrodynamic situation 
arises. The influence of frequent collisions establishes a 
quasi-equilibrium distribution which is characterized by 
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hydrodynamic variables, such as the average velocity or 
temperature which depend only smoothly on r. These 
regimes were studied in our previous publications 4 ' 8 ' 910 . 
For obvious reasons a treatment of the intermediate 
regime (I « L) presents considerable difficulties. 

In this paper we show that it is possible to study the 
intermediate regime successfully if the scattering prob- 
ability decreases substantially with each collision. The 
occurrence of a decreasing scattering probability can be 
expressed by an additional small parameter A which is 
defined as the ratio of the scattering probability of a 
nonequilibrium electron after a collision to the proba- 
bility before this event. Such situations can be found in 
a wide range of problems where non-equilibrium quasi- 
particles exhibiting some excess energy, above thermal 
equilibrium, collide with equilibrium quasi-particles and 
lose a considerable amount of this energy due to this 
event. Usually, the probability for inelastic scattering 
decreases drastically with decreasing energy of the quasi- 
particles, because the statistical weight decreases for both 
the possibility for finding a collision partner and for 
finding empty final states. Electron-electron, electron- 
phonon or phonon-phonon collisions at temperatures be- 
low the Debye temperature can be taken as an example. 

Our approach is based on the following concept: (i) a 
few collisions are treated exactly, while (ii) subsequent 
collisions are neglected, or considered perturbationally. 
We will refer to this approach as the "modified one- 
collision approximation" (MOC), which accounts for the 
fact that in most cases a restriction to an one collision 
approximation is sufficient to explain the general features 
of the scattering experiments. In the present work, a full 
account of the MOC approximation is given and the re- 
sults are applied to analyze experimental data obtained 
from electron beam experiments. It is shown, that appli- 
cability of the MOC approximation to the scattering of 
high-energy electrons substantially extends the descrip- 
tion of spectroscopic data from electron-electron colli- 
sions experiment* 4 . The non-monotonic dependence of 
the signal of scattered electrons on the beam energy can 
be explained within the MOC approximation. This ob- 
servation can be interpreted as a crossover from the bal- 
listic to multi-collision (hydrodynamic) flow regime. 



II. MODIFIED ONE-COLLISION 
APPROXIMATION 

Let us consider non-equilibrium electrons at a charac- 
teristic energy eg, measured from the Fermi level Ef, and 
larger than the temperature T of the equilibrium electron 
gas, £o 3* T. In a first stage of relaxation these highly 
non-equilibrium electrons are scattered by equilibrium 
electrons thus generating a new non-equilibrium system 
but with already much lower excess energy. This stage 
of relaxation was named the "pretemperature" stage in 
Refi This "pretemperature" stage ends when the quasi- 
particles energy becomes comparable to the bath temper- 



ature T. The most important fact for the "pretempera- 
ture" stage is the following: If the energy e < £f, then 
the probability of secondary collision is an order of mag- 
nitude lower than for the first onat Indeed, the energy e 
of a non-equilibrium electron after collision with an equi- 
librium electron will be distributed between three states 
in phase space: 

p + p'-p" + p'"; 

e= \e'\+e" + e'" . 

For T < e, we have e' < and e" , e'" > due to Pauli 
principle, i.e. the state p' is below the Fermi level while 
the final states p" and p'" are well above. The energy 
of the initial electron is redistributed between three final 
states. On average, the energy of each state after the 
first collision can be expressed as £\ ~ £o ■ A, with A 2 w 
1CP 1 (Note, the numerical value of A depends on the 
specific type of problem considered)^. Since the energy 
dependent mean- free path 1(e) cx 1/e 2 , this means that 
1(e) of non-equilibrium quasi-particles will increase by 
about an order of magnitude, which implies that l(e\) ~ 
A~ 2 • l(eo) l(eo)- This may lead to a situation where 
the probability for a secondary collision is low while the 
probability for the first collision is close to unity in a 
system of size L: l(si) ^> L > l(eo). In this case, it 
is a natural choice to consider the first collision exactly 
while the secondary collision is treated by perturbation 
theory, which leads to an extension of the validity of the 
one-collision approximation up to L ~ 10 l(eo). 

In the following, we use this approach for the case of 
propagation of nonequilibrium electrons in an electron 
beam injected into a 2DEG via a quantum point-contact. 
For steady state conditions, the operator d/dt on the left 
side of Eq. (QJ must be replaced by (vd/dr) + (Fd/dp). 
For further consideration it is convenient to use Carte- 
sian coordinates. The x-axis is along the beam direction, 
while the origin is located at the center of the electron 
beam injector (a quantum point-contact in an actual ex- 
periment). Eq. Q will be solved by integrating along 
the trajectory. We define t as the time of motion of 
the particle along a given trajectory. R(r, p, t) is the 
(Lagrange) coordinate of the particle at a time t. At 
the time of observation (taken as the origin of the time 
counting), the particle is located at the point r with a 
momentum p. The momentum of the particle at a time t 
is P(r, p, t) = mV = m(dll/dt). The particle trajectory 
R(r, p, t) is uniquely defined for given applied fields. The 
solution of the equation X(r, p, to) = defines the start 
time to of the electron injection into the system at x = 0. 
The solution of Eq. JIJ is given as 

oo 

/ = £/ (n) , 

n=0 
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/("+ 1 )(r,p)= j\t J dp'lv Pp ,-5(P-p') H/ (n) (R,P') 



dt{0-v)f n \ 



(2) 



In this expression each n-th term of the series is ob- 
tained by n subsequent actions of the collision operator 
on the initial or unpertubated function /W) which de- 
scribes the ballistic motion of the particle: 



/(°)(r,p)=/o(yo,Po) ) 
Y = Y{r,p,t ), 
Po =P(r,p,t ). 



(3) 



Here, fo(y,p) is an injected distribution where all the 
injected electrons have a characteristic energy Eq. The 
quasi-particles, which appear as the result of a collision, 
arc described by the action of the incoming part of the 
collision operator J on fo(y, p). In other words, Ofo cor- 
respond to a particle that has been scattered once, while 
00 fo describes a subsequent second scattering event, and 
so forth. Correspondingly, the typical energies of these 
states differ: f is characterized by the energy eq, for 
Ofo the characteristic energy will be E\ ~ Sq ■ A and 
for 00 fo it will be given by £2 ~ £0 • A 2 etc. Con- 
sequently we find that v~ 1 uOfg ~ I (si)Pfo and that 
Ov~ l vf$ ~ l^ 1 (eo)Ofo. Therefore, as a main approxi- 
mation in the parameters Z(eo)/i(ei), L/l{e\) <C 1, it is 
possible to neglect all terms in the series of Eq. JI1J that 
include or v to an order higher than one, where v is 
placed to the left of 0. After making this approximation, 
the summation of the remaining terms can be carried out 
easily. Thus, the main idea consists in a partial summa- 
tion of the infinite series to a first or subsequent approx- 
imation, as required in order to reach sufficient accuracy 
for an adequate description of the experiment. 

For a stationary beam in a magnetic field, a Green's 
function of the "MOC-1" (i.e. the solution of Eq. (JTJ for 
fo{y, p) = 5(y)8(p — Po), that takes into account the first 
collision exactly) has the following form: 



9(r,p) =9b + Yl ^ pp ' IvJxV'l 6XP ~ ' 

(4) 

Here gb is the distribution of the ballistic part of the 
electrons, i.e., the solution of Eq. in the absence of 
the integral operator 0. t' is the time of collision and 
to is the starting time. P' = mV and Pi = mVi 
are the momenta before and after collision, respectively, 
while Vo = Po/m. The values of t, to, R(r,p, t'), 
P' = P(0,p o ,t o - t'), Pi = P(r,p,i'), and p are de- 
termined through the equations of motion by the initial 
momentum po and location r. The equations of motion 
are reduced to R(r, p, t') = R(0, p , to—t'). The presence 
of walls in the experimental setup, see Fig. ^ excludes 
solutions with a full cyclotron orbit. 



Note that the MOC approximation is applicable in 
principle under conditions, where a pre-temperature 
stage occurs for the non-equilibrium quasi-particles. In 
other words, it holds when the following hierarchy of 
mean-free paths is given: 



l(e ) « l{ei) « l(e 2 ) « l(e 3 ) « ... < l(T), (5) 

where is the energy after the i-th collision and 1(sq) < 
L <C l(T); for L <C l(eo) the inequality corresponds to 
the quasi-ballistic regime, where perturbation theory is 
applicable. If L corresponds to the n-th segment of this 
hierarchy, i.e. l(e n ) < L < l(e n +i), it is necessary to 
construct a "MOC-n" approximation by treating the first 
n collisions exactly; which leads to relative corrections of 
the order of A 2 . The approach breaks down when the 
quasi-particle energy is comparable to T. As an example, 
the exact result for the first two collisions ("MOC-2"), for 
the spatially homogeneous case and fo oc 5{e — Sq), reads: 



/ = { exp(-v t) + (v - v) 1 (exp(-- — ^7—) 
L V 1 — V Vn 



— exp(— v t)j + v Q 
+ (>o - vy 1 (cxp(-v t) - 



v~ x (\ - exp(-i/i)) 



vt 



exp( 



))]*}fo, 



1 - v/VQ ' 

v a = u(e ). (6) 



This formula is valid for times t < 1(s2)/vf- Thus, just a 
few steps lead to a high accuracy. It turns out, that for 
the experiments discussed below, just one step is enough 
to produce satisfactory results. 



III. EXPERIMENT 

In order to investigate the electron-electron scattering 
properties in a 2DEG, experiments are needed which al- 
low to distinguish between scattered and non-scattered 
particles, and reveal the details of scattered signal. 
Therefore, we used the experimental setup of the electron 
spectroscopy method proposed in Refii (see Fig.0. The 
sample is based on a conventional Si-modulation doped 
Ga-As(Al,Ga)As heterojunction with a carrier density of 
n f» 2.8 • 10 11 cm^ 2 and an impurity mean-free path of 
li m ~ 20 /im. The scheme of the experiment is shown 
in Fig. ^ Schottky gates define quantum point-contacts 
in the 2DEG and are adjusted such that the injector and 
detector quantum point-contact (QPC) are carrying only 
one conducting mode. The distance between the injec- 
tor and detector QPC is L = 4 um. Ohmic contacts are 
indicated by crossed squares in Fig. ^ The sample was 
cooled down to T = 0.2 K. An ac-modulated dc-bias volt- 
age applied between ohmic contacts 1 and 2 was used to 
inject the electrons into the 2DEG area between injector 
and detector. The amplitude of the ac modulation 5Vi 
was much smaller than the dc component, 6Vi <C Vi. The 
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FIG. 2: The experimentally measured signal on the detector 
as a function of magnetic field for different energies of injec- 
tion: (solid) - 0.8 mV, (dash) - 1.6 mV, (dot) - 2.6 mV, 
(dash-dot) - 3.4 mV. 



ac-component of the detector signal, i.e. the voltage that 
occurs between contacts 3 and 4 (Vd = V34) is measured 
with a lock-in technique. With this technique only the 
signal of electrons with energy Vi is measured, while the 
whole e-beam consists of electrons with energies in the 
range of ... Vi- (Note that, throughout this paper the 
electron charge e is set equal to one e = 1.) 

The magnetic field dependence of the detector signal 
6Vd{B, Vi) has been studied for a wide range of energies 
Vi = . . . 9 mV. An example is given in Fig. [5] for differ- 
ent injection energies Vi. This figure display the evolu- 
tion of an mono-energetic electron beam with increasing 
injection energy Vi, which reflects the energy-dependent 
electron-electron scattering. From the half-hight width 
of low energy curve (1) (Vi = 01. mV) it is possible to ex- 
tract the opening angle of the collimated electron beam: 
<F « 12°, which is identical for injector and detector when 
the QPCs are equally adjusted to one conducting mode 
(Gqpc = 2/h). 

The contribution of scattered particles to the detector 
signal 5V d s can be easily extracted following the approach 
of Ref^: 



SVJ(B, V) = 6V d (B, V,)-exp arcsin —JsV$(B), 

(7) 

here, l(Vi) = Anhv F (e F /Vi) 2 / 'e F (lne F /Vi+\n2q TF /k F + 
1/2) is the energy dependent e-e scattering length accord- 
ing to Ref£, where <7tf is the Thomas-Fermi screening 
wave vector; r c is the cyclotron radius; v F is the Fermi- 
velocity. The second term on the right-hand side of this 
expression corresponds to gb defined in Eq. |0J, 5V^(B), 
represents the signal that would be observed in the ab- 
sence of scattering. A good approximation for this value 
can be obtained from the signal measured for the low- 
est excitation voltage SVd(B, Vi) at Vi — 0.1 mV. In this 
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FIG. 3: The maximum of scattered particles signal as a func- 
tion of injection energy. Empty square markers correspond to 
the experimental values. The thick line represents the theory. 
The dashed line represents theoretical calculations outside the 
range of applicability of the approximation used here. 

cases, L/l ~ 10 2 and electron-electron collisions can be 
ignored. Note that for injection energies Vi > 3.5 mV 
« £f/3 the above equation for l(V) yields a large inaccu- 
racy. However, this can be ignored because l(Vi) becomes 
so small that the subtraction of the ballistic contribution 
to the signal leads only to corrections of less than one 
percent. 



IV. DISCUSSION 

In Fig. |3 the maximum value of the magnetic 
field dependent signal due to the scattered electrons, 
SV^ ax (Vi) = maxs 6VJ[(B, Vi) is presented as a function 
of the injection energy (Vi). This dependence shows a 
distinct non-monotonic behavior, i.e. a rapid increase 
at small Vi, a distinct maximum at Vi ~ 2 mV and an 
abrupt decrease for larger injection energies up to Vi ~ 8 
mV, followed by a weak increase. 

The presence of a maximum and the subsequent de- 
crease in the energy dependence of SV^^iVi) can not be 
explained in the context of a simple perturbation theory 
applied to the collision integral, because in the regime 
where the maximum is observed the electron-electron 
mean-free path is already smaller than the distance be- 
tween the injector and detector QPC. Perturbation the- 
ory can only be applied to describe the increase of the 
signal for small energies Vi . The presence of a strong de- 
crease appears at first sight to be unexpected, because in 
the multi-collision regime that is entered upon increasing 
Vi, a decrease of the ac-signal from scattered particles is 
not anticipated. In this multi-collision regime the ballis- 
tic beam is already absent but the motion of the electrons 
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can be represented by a hydrodynamic drift. 

However, the behavior shown in Fig.|2|can be explained 
in the context of the modified one-collisional approxima- 
tion described above by taking into account the specific 
properties of e-e scattering in two-dimensionsi. A calcu- 
lation shows, that by taking into account the first colli- 
sion exactly, Eq. (0J leads to the following expression for 
the detector signal due to scattered particles 

8Vi oc vK{^)G{ij)) exp ( - j ^.J^ ) ,i/f = 2 arcsin — , 

(8) 

where ip is the scattering angle and G(ip) = 
X*$/2 d<pg{ip ~ f) i s the angular distribution function of 
e-e scattering, g{ip) — mv^ 1 J de v eea (ip) , averaged over 
an interval <I>; the weighting coefficient K is K m 1/$ for 
tp < $ and K l/sin-0 for $ < ip < n/2. (It is suffi- 
cient to consider only the case where the cyclotron radius 
r c is smaller than L.) It should be noted that Eq. © is 
valid for all ratios of $ and the angular width of g(ip). 
The expression Ltpj\ sin(-0/2) in the exponent of Eq. (jHJ) 
represents the length of the electron trajectory from the 
injector to the point O where the scattering takes place 
(see Fig. 0. From Eq. JSJ) it follows, that the signal 
due to scattered particles is directly proportional, on the 
one hand, to the e-e scattering probability, or the fre- 
quency of the collisions v and, on the other hand, it is 
proportional to the probability of an electron to reach 
the scattering point O ballistically. As has been shown 
earlier in^, the pre-exponential factor increases with in- 
creasing Vi at low injection energies (mainly due to the 
increase of the collision frequency v cx Vf). 

The existence of a maximum in Fig.[3]can be attributed 
to a competition between an increase of the scattering 
probability (the pre-exponential factor) and a decrease 
of ballistic transport into the neighborhood of the scat- 
tering point O (the exponent) with Vi. This explanation 
is supported by the excellent quantitative agreement be- 
tween the experimental data and a theoretical fit based 
on Eq. JHJ), where the angle ip is replaced by the experi- 
mental value ipmax , see Fig. [3] 

The good fit also justifies the use of the one-collision 
MOC approximations indicating that higher orders col- 
lisions can be neglected. For small energies, the pre- 
exponential factor determines the behaviour of the sig- 
nal due to scattered particles, which also can be described 
by standard perturbation theory applied to the collisions 
integral. Around the maximum the probability for sec- 
ondary collisions is still small and the signal is still de- 
fined by the scattering that takes place in the neighbor- 
hood of the point O - however, the probability for a first 
collision on the way between the injector and the detector 
becomes larger. This underlines the decrease of the sig- 
nal with increasing V,. As mentioned above this regime 
can not be described by ordinary perturbation theory, re- 
quing the MOC approximation. At this point we would 
like to note that the non-monotonic behavior of the ac- 



current component which has been discussed above can 
not be explained by nonlincarities due to heating effects 
which are expected to appear at approximately the same 
values of Vi in de-measurements^. It can be shown, that 
such an effect gives only a small correction of the order of 
T(Vi)/Vi to the ac-component of the signal, where T(Vi) 
is the temperature due to the heating effect. 

Finally, we may only speculate about the origin of the 
weak increase of the detector signal which is observed for 
the highest injection energies Vi. As the injection en- 
ergy approaches Sf, the potential drop across the injec- 
tor QPC becomes so large that the second subband of the 
QPC may become populated. This enhances the injector 
conductance, and, consequently, the detected signal. For 
these energies Eq. is no longer valid (see Refi^). A 
second reason could be transition to the two-collision (or 
multi-collision) regime as the mean-free path is less then 
10~ 1 L at these energies. 



V. CONCLUSION 

In summary, we have shown in this paper that it is pos- 
sible to successfully investigate the evolution of a high- 
energy non-equilibrium quasi-particles distribution with 
the help of the " modified one-collision approximation" . 
This method was used to explain novel measurements of 
the non-monotonic energy dependence of the signal of 
scattered electrons in a 2D system. We have obtained 
an excellent quantitative agreement between the experi- 
mental data and a theoretical calculation. 

In conclusion we would like to emphasize that the 
MOC approximation may be of interest for a wide range 
of problems independent of the nature of the scattering 
particles; the only restriction is that the scattering prob- 
ability should strongly decrease after one or several colli- 
sions. Along with the case of electron-electron scattering, 
this is a typical situation for the majority of inelastic 
scattering processes, such as electron-phonon, phonon- 
phonon etc. Let us discuss as an example, electron- 
phonon scattering in the case when the electron energy 
is high enough e T, but is low in comparison with the 
Debye temperature (the electron energy is measured from 
the Fermi level) . In this case a phonon with an energy of 
the order of e will be emitted as a result of the electron- 
phonon scattering event and the corresponding mean-free 
path is I cx e~ 3 . In contrast to electron-electron scatter- 
ing the energy is redistributed here between two states ac- 
cording to the momentum and energy conservation laws, 
Po = Pi + q and So = £i + w. Here po and e\ are the 
quasimomentum and energy before the collision, and af- 
ter the collision they are pi, and s\, respectively; q is 
quasimomentum of the emitted phonon and u> its energy. 
Since, on average, the energy is redistributed equally be- 
tween the scattering product lu ~ e\, we also have here 
the small parameter l(eo)/l(ei) ~ 2~ 3 ~ 0.1. The case 
of three-particles high-energy phonon-phonon scattering 
processes differs only by the energy dependence of the 
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mean-free path. 

Furthermore, the inelasticity, on which the approach 
developed in this paper is based, is likely to be relevant 
not only to the energy redistribution but also to some 
other quantities; this includes, for example, the spin in 
cases where the probability of scattering decreases signifi- 
cantly after a spin- flip scattering event. This case may be 
realized for electron scattering from spin-polarized impu- 
rities (there are no spin subbands in this case) and if the 
exchange interaction depends significantly on the mutual 
orientation of the electron spin and the impurity spin. 
Let us neglect the interaction between electrons of one of 
the spin directions (spin-up, for example) and the impu- 
rities. Then the spin relaxation process is completed just 
after the first spin-flip scattering of the spin-down elec- 
trons. Thus, it is straightforward to describe this process 
in the framework of the modified one-collision approxi- 
mation. 

The method we proposed in this work can be very ef- 
fective also for point-contact spectroscop y 11 ' 12 . The rea- 
son is the following. The current- voltage characteristic of 
the electric circuit containing a microconstriction is de- 
termined by the back-scattering due to inelastic electron 
scattering processes just in the region of the microcon- 
striction. If 1(e) < L (L is the length of the constriction) 



the current-voltage characteristic is determined only by 
one, or few, first collisions, the probability of which is not 
small and, thus, it can not be taken into account using 
a perturbation theory (as it was done in the study de- 
scribed in Ref. 12 for the case when 1(e) S> L). The same 
also holds for the nonlinear thermoelectric phenomenon 
in point-contacts^. The influence of the first collisions 
on the point-contact spectrum may be calculated exactly 
using equations of the general form of Eq. (6) because in 
point-contact spectroscopy the excess energy of the elec- 
tron is high, £>T. 
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